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A new method of pseudo-state discretization is proposed for the method of continuum discretized 
coupled channels (CDCC) to deal with three-body breakup processes. We propose real- and complex- 
range Gaussian bases for the pseudo-state wave functions, and show that they form in good ap- 
proximation a complete set in the configuration space which is important for breakup processes. 
Continuous S-matrix elements are derived with the approximate completeness from discrete ones 
calculated by CDCC. Accuracy of the method is tested quantitatively for two realistic examples, 
d+ 58 Ni scattering at 80 MeV and 6 Li+ 40 Ca scattering at 156 MeV with the satisfactory results. 
Possibility of application of the method to four-body breakup processes is also discussed. 

PACS numbers: 24.10.Eq, 25.45.De, 25.60.Gc, 25.70.Ef 



I. INTRODUCTION 

The method of continuum discretized coupled channels 
(CDCC) has been successful in describing nuclear reac- 
tions inclu ding weakly bound projectiles [E IE IE El IE IE 
HE IE EE ETUI- CDCC has been attracting much at- 
tention since the advent of experiments with radioactive 
beams, because projectile breakup processes are essential 
to many of such reactions. CDCC plays an important role 
in the spectroscopic studies of radioactive nuclei through 
the nuclear reactions involving such nuclei. 

In CDCC for reactions with a projectile consisting of 
two fragments, the states of the projectile are classified 
by the linear and the angular momenta, k and i, of rela- 
tive motion of the two fragments of the projectile, which 
are truncated by k < fc max and I < f max . The trun- 
cation is the most basic assumption in CDCC, and it 
is confirmed that calculated S-matrix elements converge 
for sufficiently large &; max and £ max 0, IE E3 • It has been 
shown that CDCC is the first-order approximation to the 
distorted Faddeev equations, and corrections to the con- 
verged CDCC solution are negligible within the region of 
space in which the reaction takes place [Til ]. 

As a consequence of the truncation, the integral equa- 
tion form of the equations of coupled channels, derived 
from the three-body Schrodinger equation, has a com- 
pact kernel, satisfying a necessary condition for iterative 
solutions [14J. In practice, however, the coupled chan- 
nels equations thus obtained are impossible to be solved 
because of the continuously infinite number of coupled 
breakup channels. The problem is solved by discretizing 
the k continuum. The discretization leads the coupled 
equations to a set of differential equations with a finite 
number of channels. 
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As for the discretization, three methods have been pro- 
posed so far: the average (Av) 0,0, IE El, the midpoint 
(Mid) 0, El, and the pseudo-state (PS)[lllEII3 meth- 
ods. In the Av and Mid methods, the /c-continuum is 
divided into a finite number of bins. In the former the 
continuum channels within each bin are averaged into a 
single channel, while in the latter they are represented 
by the channel at a midpoint of the bin. It has been 
confirmed that calculated S'-matrix elements converge 
as the width A of the bins is decreased, and also that 
the two methods yield the same converged S'-matrix el- 
ements 0, |E El • From practical point of view, the Av 
method is more convenient than the Mid one, since the 
former requires less numerical works than the latter. The 
Av method therefore is most widely used in practice. 

In the PS method, the breakup states are described by 
superpositions of L 2 -type basis functions. The wave func- 
tions of such pseudo breakup state have wrong asymp- 
totic forms. For this reason, the PS method was mainly 
used in the past to describe virtual breakup processes in 
the intermediate stage of elastic scattering and (d, p) re- 
actions. Very recently, the applicability of the PS method 
for real breakup reactions for a specific set of basis func- 
tions. It was shown that the PS method well simulates 
the angular distribution of the total breakup cross sec- 
tion calculated with the Av method 0]- In the analysis, 
discrete breakup cross sections calculated with the PS 
method are transformed into continuous one, by assum- 
ing a specific form for the k distribution of the continuous 
cross section, that is, a histogram with widths estimated 
in a reasonable way. 

The purpose of the present paper is to propose a new 
PS method for projectile breakup reactions. If the basis 
functions form in the good approximation a complete set 
in the configuration space which is important for breakup 
processes, we can formulate a method of interpolation for 
generating continuous breakup S-matrix elements from 
the discrete ones provided by the PS method, following 
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the formulation of 18] for the Av method in which the 
completeness of discretized states is well satisfied if the 
discretization is enough precise. The method does not 
assume any form a priori for k distributions of continu- 
ous breakup S'-matrix elements. The method is indepen- 
dent of the type of the basis functions used for the PS 
method. Only condition is that the basis functions con- 
stitute an approximate complete set in the wide range of 
k and its conjugate coordinate r which are important for 
the breakup processes. As basis functions satisfying this 
condition, we propose two types of them. One is ordinary 
Gaussian functions |19|. which we refer to as real-range 
Gaussian function in the present paper. The other is a 
natural extension of that, complex-range Gaussian func- 
tion |2Cj . i.e, Gaussian functions with the complex range 
parameters the precise definition of which is given in a 
later section. 

In Section II, we recapitulate CDCC based on both 
the Av and the PS methods of discretization. In Section 
III, we present a method of interpolation to get contin- 
uous breakup S'-matrix elements from discrete ones cal- 
culated with the PS method, and introduce the real and 
the complex-range Gaussian basis. In Section IV, the 
validity of the present PS method is tested and justified 
for two realistic cases, d+ 58 Ni scattering at 80 MeV and 
6 Li+ 40 Ca scattering at 156 MeV. In Section V, the poten- 
tiality of the present PS method for four-body breakup 
reactions is discussed. Section VI gives a summary. 



II. THE METHOD OF CONTINUUM 
DISCRETIZED COUPLED CHANNELS 

We consider a reaction of a weakly bound projectile 
(B) impinging on a target nucleus (A). We treat a simple 
system shown in Fig. 1 in which the projectile, B, is com- 
posed of two particles (b and c) and the target is inert. 
The three-body system is described by a model Hamilto- 
nian H = JJt, c + Kr + U, where JJbc = K r + Vb c (r) and 
U = C/bA^bA) + l^cA^ca)- Vector r is the relative coor- 
dinate between b and c, R the one between the center-of- 
mass of the b-c pair and the center-of-mass of A, and rxy 
denotes the relative coordinate between two particles X 
and Y. Operators K r and Kr are kinetic energies associ- 
ated with r and R, respectively. Vb c (r) is the interaction 
between b and c. The interaction U^a {U c a) between b 
(c) and A is taken to be the optical potential for b+A 
(c+A) scattering. For simplicity, the spin dependence of 
the interactions is neglected. Furthermore, the Coulomb 
part of U is treated approximately as a function only of 
R, i.e., we neglect Coulomb breakup processes, and focus 
our attention on nuclear breakup. 

In CDCC, the three-body wave function \I/jm, with 
the total angular momentum J and its projection M on 
z-axis, is expanded in terms of the orthonormal set of 




FIG. 1: Illustration of a three-body (A+b+c) system. The 
symbol B=b+c stands for the projectile, and A is the target. 



eigenstates $ of Hbc- 

* JM (r,R) = Y,yjMMr)xtoLj{Po,R)/R 

L 

poo 

+ J2 y JM $e(k,r)xeLj(P,R)/Rdk, (1) 

l,L J ° 

where 

yjM = \i e Y e (Vr) ® i L Y L (n R )] JM - (2) 

For simplicity, we assume that the b+c system has one 
bound state ^a(r) with angular momentum £ and con- 
tinuum states <$>t(k, r) with linear momentum k and an- 
gular momentum t, both ranging from zero to infinity. 
The <$>i are real functions normalized to the <5-function 
in k 0. The projectile B is initially in the bound state. 
The coefficient xilj (X£ lj) of the expansion describes a 
center-of-mass motion of the b-c pair relative to A in the 
state with the linear and orbital angular momenta P 
(Jo) an d L, respectively. 

In CDCC, the sum over £ is truncated by I < ^ max 
and the k integral by k < /c max . For each £, furthermore, 
the continuum states from k — to fc max are discretized 
into a finite number of states, with the wave function 
of the i-th state corresponding to the momentum 
ki. Details of the discretization methods are described in 
the next section. 

After the truncation and the discretization, \&jm is 
reduced to an approximate one, 

*5m CC = E yjM®o(r)x^(Po, R)/R 

L 

! m „ N 

+ EEE^M d ^( r )^(A,i?)/i?, (3) 

1=0 i=l L 

where 

X~ t0 {Pa,R) = x~h,(Po,R), To = (0,4, J, J), 
Xj(Pi,R) - W^{Pi,R), j = (i,l,L,J). 

On the right hand side of Eq. @, the first term repre- 
sents the elastic channel denoted by 70 and the second 
one corresponds to the discretized breakup channels, each 
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denoted by 7. The weight factor W 1 depends on the dis- 
cretization method used. Each pair of momenta, (fcj, Pi), 
satisfies the total energy conservation: 

E = h 2 P 2 /2fi AB + eo = h 2 P 2 /2fi AB + et, (4) 

where eo and e< = ti 2 k 2 /2fi\ lc are the energies of the 
ground and the continuum states, respectively. Insert- 
ing Eq. J2J into the three-body Schrodinger equation, 
(H — E)^fjM = 0, leads to a set of coupled differential 
equations for x 10 (P , R) and X 7 (-Pj, R): 



d 2 
dR 2 



P 



L(L + 1) 2/iAB 



R? 



Ir 



Vry(R) 



X;^^..«/?-.v,(7V.7?) 



7V7 



h 2 



(5) 



for all 7 including 70, where ko = \/ — 2^ c e^/h and Po 
P - The coupling potentials V '(R) are obtained as 



V W {R) = (Yjt t *it{r)\U\y j £ $ iY MW (6) 

The coupled equations are solved, under the asymptotic 
boundary condition 



x^(Pi,R) 



u ( L ) (P i ,R)6 ino 




u i+) 

7,7o a L 



(Pi,R). 



Here u { L \p u R) and v£'{Pi,R) are incoming and out- for k belonging to the i-th bin, i.e., < k < fcj. 



,(+) 



(7) 



then the weight factor W 7 is given by 



W 2 = 



[Mk)] 2 dk. 



fci- 



For a bin far from a resonance, it is natural to set 



fu(k) — 1, so that W 7 = 



since <£v(fc,r) changes 



smoothly with k. On the other hand, $>e(k,r) changes 
rapidly across the resonance. One way of coping with 
this situation is to take much smaller than the width 
of the resonance so that $>t(k, r) does not change much 
within individual bins. This, however, make the num- 
ber of bins large. Alternatively, one can take a single 
bin which contains the whole resonance peak and use a 
weight function of Breit-Wigner type P, IE IE IS EH > 



fu(k) 



iT/2 



e(k) 



iT/2 



(9) 



where e(k) is a continuous intrinsic energy of the b+c 
system. The discretized intrinsic energy, ej = h 2 k 2 /2fib c , 
corresponding to each bin is obtained as k 2 = (hi + 
ki-1) 2 /A + for a non-resonance bin and e, = e res 

for a resonance one. Comparing the approximate form 
@ with the exact one in the asymptotic region 
R — > 00, it is natural to assume 



fu(k) (forAv) 



(10) 



going Coulomb wave functions with the momentum P iy 
and S^~ is the 5"-matrix element for the transition from 
the initial channel 70 to 7. 



B. The pseudo-state method 



III. DISCRETIZATION OF k CONTINUUM 

Among the three methods of discretization of the k 
continuum, the relation between the average (Av) and 
the midpoint (Mid) methods has already been clari- 
fied ^El- The present discussion therefore is focused on 
the relation between the Av and the pseudo-state (PS) 
methods. 



A. The average method 

In the Av method, the fc-continuum [0, fc max ]) for each 
£, is divided into a finite number of bins, each with 
a width An — ki — fcj-i, and the continuum breakup 
states in the i-th bin are averaged with a weight function 
fu(k) 0,0- The resultant orthonormal state is described 
as 

= 7^r [ Mk,r)f a (k)dk (for Av) , (8) 

W 1 Jki-i 



In the PS method, H\, c is diagonalized in a space 
spanned by a finite number of L 2 type basis functions. 
The resultant eigenstates can well reproduce both bound 
and continuous states within a finite region of k and 
r 15]. The k continuum is automatically discretized 
by identifying the eigenstates of positive energies with 
$i^(r). The weight factor W 1 is unity if the resultant 
discretized states <&u(r) are orthonormalizcd. Among the 
eigenstates, only low-lying states belonging to the region 
< e < h 2 k' 2 naK /2fibc are taken as breakup channels in 
CDCC equation JSJ, where fc max is the maximum k in 
the Av method. 

The CDCC equations iJJjJ thus obtained yield discrete 
breakup S'-matrix elements. If the basis functions form 
a complete set with good accuracy in the region of r and 
k which is important for breakup processes, an accurate 
transformation from the approximate breakup S'-matrix 
elements to the continuous ( "exact" ) ones is possible, as 
shown in [l^ for the Av method. The exact breakup 
T-matrix element is given by 

T^\k) = (<Mfc, r)j L (PR)y e jM\U\^j M ) ■ (11) 
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Inserting the approximate complete set {&u(r)} between 
the bra vector and the operator U in Eq. Ijlll) . and re- 
placing the ket vector by the CDCC wave function J2J), 
one obtains the following approximate relation, 



where 



and 



7.70' 



fr e s (k) = (^(fc,r)|#«(r)) 



,CDCC\ 

;m / 



(12) 



(13) 



(14) 



The last form of Eq. H12fl has been derived by replacing 
P by Pi in the spherical Bcsscl function ji,(PR), which 
is valid since the k distribution of f^ s (k) is sharply lo- 
calized at k = ki. T 7 . 7o is a CDCC breakup T-matrix 
element calculated with CDCC. Since the T 7l7o are pro- 
portional to the corresponding S'-matrix elements S. 



7,70 i 



(15) 



The "interpolation formula" (|15fl for the PS method 
agrees with the corresponding one in 01 f° r the Av 
method. Thus, the interpolation formula (|15[) can be 
used for any method of discretization, if the discretized 
wave functions constitute an approximate complete set. 
The interpolation formula l|15|) is also independent of the 
type of the basis function taken, as obvious from the 
derivation, but it is necessary that the basis functions 
form a complete set in good approximation in the region 
of r and k which is important for the breakup T-matrix 
elements. As such basis functions, we here propose two 
types; one is the conventional real-range Gaussian func- 
tions 



r'exp[-(r/oj) 2 ], 



(16) 



where a,j (J = l —n) are assumed to increase in a geomet- 
ric progression [lfij : 



a ] =a 1 (a n /a 1 )^- i y^- 1 \ 



(17) 



The other is an extension of Eq. (|16|) 20] : pairs of func- 
tions 

r^exp — (r/aj) 2 cos (b (r/aj) 2 ^j , (18a) 



r exp 



(r/djY 



{b(r/a 3 y 



which can be also expressed as (4>u + <fr* e )/2 and 
4>it)/{2i), respectively, with 

(j>it{r) = /expf-^r 2 ], rjj = (1 + ib)/a 2 , 



(19) 



(18b) 



i.e., a Gaussian function with a complex range parame- 
ter. We refer to the new basis (I18JI as the complex-range 
Gaussian basis. In Eqs. I|18l) and i|19f) aj are the same 
as in Eq. (|17f) : b is a free parameter, in principle, but 
numerical test show that b — tt/2 is recommendable. It 
should be noted that the total number of basis functions 
is 2n. Note also that the complex-range Gaussian basis 
agrees with the real-range Gaussian basis when 6 = 0. 

The complex-range Gaussian basis functions are oscil- 
lating with r. They are therefore expected to simulate 
the oscillating pattern of the continuous breakup state 
wave functions better than the real-range Gaussian basis 
functions do. Actually, the latter reproduces the contin- 
uous state <f>i(k,r) in the region < kr < 20 0], while 
the former does in the even larger region, < kr < 35. 

If necessary, one can calculate, with the complex-range 
Gaussian basis as well as the real-range Gaussian basis, 
all the coupling potentials analytically by expanding the 
potential U into a finite number of Gaussian functions. 
As far as three-body breakup reactions are concerned, 
however the analytic forms of the coupling potentials are 
not necessarily, since they can be easily obtained with 
numerical integration. Direct numerical calculation is 
difficult in the case of four-body breakup, and the an- 
alytic forms of the potentials are very useful in that case, 
as shown in section Ivl 



IV. NUMERICAL TEST OF THE 
PSEUDO-STATE METHOD 

In previous CDCC analyses 0,03, 

calculated elastic 

and breakup S-matrix elements were found to converge, 
for sufficiently large model space. In this section, we test 
the PS method by comparing the calculated S-matrix 
elements with those obtained with the Av method that 
converged within the error of 1% and hence forth called 
"exact" S'-matrix elements. The test is made for two sys- 
tems, d+ 58 Ni scattering at 80 MeV and 6 Li+ 40 Ca scat- 
tering at 156 MeV. 



A. d+ 58 Ni scattering at 80 MeV 

The model space taken in the present CDCC calcula- 
tions is I = 0,2 and fc max = 1.3 frn - . It should be noted 
that the p-wave (I = 1) breakup is negligible, because 
couplings V 11 ' (R) between odd and even parity breakup 
states contain a small factor, U = U p a(^ p a) ~ C^nA(r p A), 
in Eq. The choice of the fc max is discussed below. Ta- 
ble I shows the parameters of the potentials used; the in- 
teraction between a nucleon and the target is the nucleon- 
nucleus optical potential of Becchetti and Greenlees [2l| 
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TABLE I: Parameters of the optical potentials for n + Ni and p + Ni at the half the deuteron incident energy. We followed 
the same notation as in Ref. 21j 

system V (MeV) r (fm) a (fm) W (MeV) r w (fm) a w (fm) Wn (MeV) r WD (fin) owd (fin) 

p + 5 «Ni 44.921 1.17 0.750 6.10 1.32 0.534 2.214 1.32 0.534 

n + 58 Ni 42.672 1.17 0.750 7.24 1.26 0.580 2.586 1.26 0.580 



at half the deuteron incident energy. The interaction 
between proton and neutron is a one-range Gaussian po- 
tential, v np = vo exp[— (r/r ) 2 ] with v Q — —72.15 MeV 
and ro = 1.484 fm, which reproduces the radius and the 
binding energy of deuteron. 

In the Av method the weight function is taken as 
fuik) = 1, since the projectile (deuteron) has no res- 
onance state. The model space that gives convergence 
within error of 1% turns out to be An = 1.3/30 fm -1 , 
as mentioned above, and fc max = 1.3 fm -1 . The resulting 
values of fc max and Aj^ are different from those used in the 
previous analysis [l3j ; the main purpose of Ref. [13J was 
to show that the convergence of the CDCC solution was 
obtained within a model space of practical use and that 
the converged solution satisfied an appropriate bound- 
ary condition. The model space taken there, £; max = 1.0 
fm -1 and An — 1/8 fm -1 , are indeed enough for the 
elastic S-matrix elements and the dominant part of the 
breakup ones with the smaller k, therefore the elastic 
cross sections and the total breakup cross sections are 
well reproduced. However, the model space is found to 
be insufficient to obtain the "exact" S'-matrix elements 
in the high k region around 1.0 fm -1 , hence we take 
here fc max = 1.3 fm -1 and An = 1.3/30 fm -1 mentioned 
above. 

In the real-range Gaussian PS method, a similar con- 
vergence is found, when the number of breakup chan- 
nels, iVps, is 18 for both s- and d- waves. The number 
is even smaller when the complex-range Gaussian ba- 
sis is taken: Nps is 16 for s-wave and 17 for d-wave. 
The basis functions finally obtained have parameter sets 
(ai = 1.0, a n — 30.0, n = 30) for real-range Gaussian 
basis and (a a = 1.0, a n = 20.0, 2n = 40,6 = tt/2) for 
complex-range one. For both of them Npg is smaller than 
the number of basis functions. High- lying states with 
k > fcmax, which are obtained by diagonalizing H\> Cl do 
not affect the breakup S'-matrix elements with k < fc max , 
because the coupling potentials between the two k regions 
are weak. 

Figure 2 shows the discrete momenta hy translated 
from the eigenenergies ea for the real- and complex-range 
Gaussian bases. One sees that for the real-range Gaus- 
sian basis, the discrete momenta are dense in the smaller 
k region and sparse in the larger k one. This distribu- 



(a) ■ 
1.3 



(b) 



-2.22 MeV 

real-range Gaussian basis 



-2.22 MeV 

complex-range Gaussian basis 



FIG. 2: Discretized momenta for real-range (a) and complex- 
range (b) Gaussian bases for deuteron. In each panel, the left 
(right) side corresponds to the s-state (d-state). The hori- 
zontal dotted line represents the maximum momentum fc max 
taken to be 1.3 fm -1 . 
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FIG. 3: Angular distribution of the elastic differential cross 
section as a ratio to the Rutherford one for d + 58 Ni scatter- 
ing at 80 MeV. Results with the Av method, the real- and 
complex-range Gaussian PS methods are represented by the 
solid, dashed and dotted lines, respectively. 
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FIG. 4: The squared moduli of breakup S'-matrix elements as a function of k at the grazing total angular momentum J — 17 
for d + 58 Ni scattering at 80 MeV. The panel (a), (b), (c) and (d) shows the result for the s-state with, and the d-states with 
L — 15, 17 and 19 respectively. In each panel, the dashed (solid) line represents the result of the real-range (complex-range) 
Gaussian PS method. The step line is the result for the "exact" S'-matrix element calculated the Av method. 



tion is not so effective in simulating the k continuum, in 
the higher k region in particular. For the complex-range 
Gaussian basis, on the other hand, the discrete momenta 
are distributed almost evenly. A similar sequence of the 
k\ is also seen for the case of the transformed harmonic 
oscillator basis of Ref. 0, ^ . Such a sequence of ki is 
close to that in the Av method. Thus, the complex-range 
Gaussian basis, as well as the transformed harmonic os- 
cillator, is well suited for simulating the k continuum in 
the entire region < k < k max . 

For the elastic S-matrix elements, both the real- and 
complex-range Gaussian PS methods well reproduce the 
"exact" one calculated with the Av method, as confirmed 
in Fig. 3 for the differential cross section. The three types 
of calculations, the real-range Gaussian PS, the complex- 
range Gaussian PS and the Av methods, yield an iden- 



tical cross section at all scattering angles. Thus, both of 
the PS methods proposed here are useful for treating the 
breakup effects on the elastic scattering. 

Figure 4 shows the result for breakup S'-matrix ele- 
ments at the grazing total angular momentum J = 17, 
as a function of k. The real-range Gaussian PS method 
(dashed line) well simulates the exact solution calculated 
with the Av method (step line) in the lower k region that 
corresponds to the main components of the breakup S'- 
matrix elements, but inaccurate in the higher k region 
around k = 0.8 fm" 1 . The deviation at higher k stems 
from the fact that real-range Gaussian basis poorly re- 
produces the continuum breakup state <&e(k,r) at the 
higher k. Figure 4 shows that this problem can be solved 
by using the complex-range Gaussian basis (solid line) 
instead. 
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(a) 



i 1.0 



(b) 



s-state d-state 



1 .47 MeV 

real-range Gaussian basis 



s-state d-state 



1 .47 MeV 
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FIG. 5: The same as in Fig. 3 but for 6 Li; fc max is taken to 
be 2.0 fm -1 . The horizontal dashed line corresponds to the 
border momentum between the resonance and non-resonance 
parts used in the Av method. 
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FIG. 6: Angular distribution of the elastic differential cross 
section (Rutherford ratio) for 6 Li + 40 Ca scattering at 156 
MeV. The results of the complex-range Gaussian PS method 
and the approximate treatment of the resonance of 6 Li, i.e., 
the conventional Av method with the weight factor of Breit- 
Wigner type, are shown by the dashed and dash-dotted lines, 
respectively. The solid line is the "exact" solution calculated 
by the Av method with dense bins and the dotted line is the 
result of Watanabe model, i.e., without breakup effects. 



B. 6 Li+ 40 Ca scattering at 156 MeV 

Characteristic to this scattering, the projectile ( 6 Li) 
has d-wave triplet resonance states (3 + ,2 + ,l + ). For 
simplicity, we neglect the intrinsic spin of 6 Li, following 
Refs. pHEEl- Then the projectile has only one d-wave 



resonance state with e res = 2.96 MeV and T = 0.62 MeV. 
Obviously the energy and the width do not reproduce ex- 
perimental data, but at least the elastic cross section of 
6 Li is not affected much by the neglection of the spin [2j| . 

In this scattering, the three-body system consists of 
deuteron, a and 40 Ca. The interactions between each 
pair of the constituents are the optical potential of 
a+ 40 Ca scatterin g at 104 MeV [24|, that of d+ 40 Ca scat- 
tering at 56 MeV (25|, and v a d = voexp [— (r/r ) 2 ] with 
v = -74.19 MeV and r = 2.236 fm. Table II shows the 
parameters of the optical potentials. 

The model space sufficient for describing breakup pro- 
cesses in this scattering is fc max = 2.0 fm -1 and £ max = 2; 
the model space is composed of two fc-continua for £ — 
and 2. Since there exists a resonance in £ = 2, the d- 
wave fc-continuum is further divided in the Av method 
into the resonant part [0 < k < 0.55] and the non- 
resonant part [0.55 < k < 2.0]. The k continuum of 
&i,£=2(k, r) in the resonant part varies rapidly with k. 
The Av method can simulate the rapid change by taking 
fi.e=2(k) — 1 with bins of an extremely small width. In 
fact clear convergence is found for both the elastic and 
the breakup 5*-matrix elements, when the resonance part 
is described by 30 bins and the non-resonance part of 
the d-wave and the s-wave fc-continua by 20 bins. An- 
other Av discretization, which has been widely used as 
a convenient prescription 0, IE IE IE EJ- * s a l so ma de 
for comparison, in which the resonance region is repre- 
sented by a single state with the weight factor of Breit- 
Wigner type jnjl. The two sorts of Av discretization 
are compared with the real- and complex-range Gaus- 
sian PS methods. With the PS methods, convergence of 
the jS-matrix elements is found with 21 s-wave breakup 
channels and 22 d-wave ones. The level sequence of 
the resulting discrete eigenstates is shown in Fig. 5 for 
both the basis functions. The level sequences have the 
same properties as in Fig. 2. The parameter sets of 
the basis functions, finally taken in the PS methods, are 
(ai = 1.0, a n = 20.0, 2n = 40, b = tt/2) for the complex- 
range Gaussian basis and {a\ — 1.0, a n — 30.0, n = 30) 
for the real-range one. 

Figure 6 shows the differential cross section of the elas- 
tic scattering. The result with the precise Av discretiza- 
tion based on dense bins, considered to be the "exact" 
solution, is denoted by the solid line. The dotted line 
represents the result of the Watanabe model, i.e., with 
no breakup channels. The conventional Av discretiza- 
tion, based on the weight factor of Breit-Wigner type 
(dash-dotted line), well describes the breakup effect, par- 
ticularly at very forward angles (0 < 20°), but deviates 
considerably from the "exact" solution at larger angles 
(9 > 30°). The complex-range Gaussian PS discretiza- 
tion (dashed line) well reproduces the "exact" solution 
with a number of channels being suitable for practical 
use. The real-range Gaussian PS method gives just the 
same result as the complex-range one. 

Figure 7 represents breakup ^-matrix elements at graz- 
ing total angular momentum J — 43. The real- and 
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TABLE II: The same as in Table I but for a + 40 Ca at 104 MeV and d + j8 Ca at 56 MeV. 
system Vb (MeV) r (fm) a (fm) Wo (MeV) r w (fm) a w (fm) Wo (MeV) r W D (fm) a W D (fm) 
a + 40 Ca 219.30 1.21 0.713 98.8 1.40 0.544 

d + 40 Ca 75.470 1.20 0.769 2.452 1.32 0.783 9.775 1.32 0.783 
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FIG. 7: The same as in Fig. 4 but for 6 Li + 40 Ca scattering at 156 MeV. The corresponding grazing total angular momentum 
is 43. The step line is the result of the Av method with dense bins. Note that the difference between the results of the real- 
and complex-range Gaussian PS methods is not visible since it is less than about 1%. 



complex-range Gaussian PS discretization well reproduce 
the "exact" solution calculated by the Av discretization 
with dense bins. The results of the two PS methods turn 
out to coincide within the thickness of the line. The res- 
onance peak can be expressed by only 8 (12) breakup 



channels in the complex-range (real-range) Gaussian PS 
method, while the corresponding number of breakup 
channels is 30 in the Av method, as mentioned above. 
Thus, one can conclude that the real- and complex-range 
Gaussian PS methods are very useful for describing not 
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only non-resonant states but also resonant ones. 

V. DISCUSSIONS ON FOUR-BODY BREAKUP 
REACTION 

In the past CDCC calculations the projectile was as- 
sumed to be a two-body system, dealing only with three- 
body breakup reactions. In this section, we investigate 
the applicability of CDCC to four-body breakup reac- 
tions of the projectile consisting of three particles, b+c+x 
(Fig. 8). The Av method needs the exact three-body 
wave functions being impossible to obtain. We can cir- 
cumvent this problem with the present PS method; one 
can prepare an approximate complete set {$^} by di- 
agonalizing the Hamiltonian of the projectile in a space 
spanned by a set of basis functions of L 2 type. With 
{<&u } as the wave functions of the breakup channels, one 
can obtain an approximate total wave function v]/ CDCC by 
solving CDCC equations 101 . Inserting i|f CDCC into the 
exact form of breakup T-matrix elements in place of the 
exact total wave function, one reaches an approximate 
form: 

T 4 = <e l ( p - R + k - r +^)|[/ 4 |* CDCC ) R , r , y , (20) 

where U4 is the sum of all interactions in the four-body 
system (A+b+c+x) , r and y are and k and q are, respec- 
tively, two Jacobi coordinates of the three-body (b+c+x) 
system and momenta conjugate to them. The accuracy of 
Eq. (|2U[1 depends on how complete the set $^ is within 
the region (0 < r < r max ,0 < y < y ma x) which is im- 
portant for the breakup process considered. An impor- 
tant advantage of the use of the real- and complex-range 
Gaussian bases is, as mentioned before, that analytic in- 
tegrations over r and y can be done in Eq. (|20|l . by ex- 
panding U4 in terms of a finite number of Gaussian basis 
functions. This makes the derivation of T4 feasible. Anal- 
yses based on this formulation are of much interest as a 
future work. 



X 




breakup processes, in which the discretization of the k 
continuum is essential. In this paper, we propose the 
new method of pseudo-state (PS) discretization which 
can be used not only for virtual breakup processes in 
elastic scattering but also for breakup reactions. First 
we show that an accurate transformation from the dis- 
crete breakup S-matrix elements calculated with the PS 
method to continuous ones is possible, since the PS basis 
functions can form in the good approximation a complete 
set in the finite region of r and k which is important for 
the breakup processes. As bases satisfying the complete- 
ness, we propose the real- and complex-range Gaussian 
bases. Both bases can treat virtual breakup processes in 
the elastic scattering with high accuracy, i.e., with the er- 
ror of calculated cross sections less than 1%. For breakup 
processes, the complex-range Gaussian basis is accurate 
throughout the entire region of the k continuum. The 
real-range Gaussian basis also keeps a good accuracy for 
the dominant part of breakup 5*-matrix elements with 
the lower k, although it is partially inaccurate for the 
higher k region. Thus, both bases can be used for real- 
istic analyses of elastic scattering and projectile breakup 
reactions. 

The present new PS method has at least two ad- 
vantages over the widely used momentum bin average 
method. One is that it does not need the exact wave 
function of the projectile over the entire region of r. This 
is important from a theoretical point of view. The other 
advantage is that with the real- and complex-range Gaus- 
sian bases one can calculate all the coupling potentials 
semi-analytically, which is very useful in actual calcula- 
tions. Furthermore, if the projectile has resonances in its 
excitation spectrum, the new method discretizes the com- 
plicated spectrum with a reasonable number of the ba- 
sis functions, without distinguishing the resonance states 
from non-resonant continuous states. These advantages 
of the new method are extremely helpful, sometimes even 
essential, in applying CDCC to four-body breakup ef- 
fects of unstable nuclei such as 6 He and 11 Li. Actually, 
a CDCC analysis of four-body breakup effect on the 6 He 
elastic scattering is in progress, and the result of the anal- 
ysis will appear in a forthcoming paper. 



FIG. 8: Illustration of a four-body (A+b+c+x) system. The 
projectile consists of b, c and x, and A is the target. 
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